Quantum Coulomb glass — Hartree-Fock approximation versus exact diagonalization 
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We investigate the behavior of disordered interacting electrons in the insulating regime. Our study 
is based on the quantum Coulomb glass model which is obtained from the classical Coulomb glass by 
adding hopping matrix elements between neighboring sites. We use two different numerical methods, 
viz. a Hartree-Fock approximation and an exact diagonalization and compare the results for the 
tunneling density of states and the localization properties in order to determine the range of validity 
of the Hartree-Fock method. We find that the Hartree-Fock method gives a good approximation 
for the density of states for all energies but represents the localization properties correctly close to 
the Fermi level only. Some consequences for the localization of disordered interacting electrons are 
discussed. 
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The physics of disordered interacting electrons has 
been a subject of great interest within the last two 
decades. Most work, both experimental and theoreti- 
cal has concentrated on the metallic regime where ex- 
periments are easier to carry out and theoretical studies 
can be based on estabHshed perturbative methods [Q. 
In comparison, the insulating regime has seen much less 
activities. Experimentally, this is due to the fact that 
transport properties vanish in an insulator at zero tem- 
perature. On the theoretical side, the main reason is that 
perturbative methods cannot be applied since the insu- 
lating (classical) limit |^,|| itself is a complicated many- 
particle problem. The prototype model in the classical 
insulating regime is the Coulomb glass model Q which 
describes the electrons as classical point charges. 

Recently, the attention has focused on the quantum 
insulating regime closer to the metal-insulator transition 
(MIT) where the description in terms of classical point 
charges becomes questionable. These new studies try to 
address two main problems: (i) one wants to gain an un- 
derstanding of the MIT itself by approaching it from the 
insulating side which is complementary to the usual ap- 
proach based on perturbation theory around the metallic 
system; (ii) one wants to know whether the key proper- 
ties of the insulating limit as, e.g., the Coulomb gap in 
the single-particle density of states (DOS) carry over 
to the regime close to the MIT where most of the exper- 
iments on insulators are performed. 

The prototype model for the quantum insulating 
regime is the quantum Coulomb glass model It is 

defined on a regular hypercubic lattice with N = L'^ {d is 
the spatial dimensionality) sites occupied by KN (spin- 
less) electrons (0 < if < 1). To ensure charge neutrality 
each lattice site carries a compensating positive charge of 
Ke. The Hamiltonian of the quantum Coulomb glass is 
obtained from the classical Coulomb glass by adding hop- 
ping terms of strength t between nearest neighbor sites. 



It reads 



+ - + \ ^(n. - K){nj - K)U,j (1) 

where c\ and q are the electron creation and annihilation 
operators at site z, respectively, and (ij) denotes all pairs 
of nearest neighbor sites, rii is the occupation number of 
site i and ji is the chemical potential. The Coulomb inter- 
action Uij = jvij remains long-ranged since screening 
breaks down in the insulating phase. The random po- 
tential values ipi are chosen independently from a box 
distribution of width 2Wo and zero mean. 

In a previous paper |^] we have investigated the quan- 
tum Coulomb glass by means of a Hartree-Fock (HF) 
approximation. Within this method the interaction was 
treated at HF level and the arising self-consistent disor- 
dered single-particle problem was diagonalized numeri- 
cally. This method enabled us to study comparatively 
large systems of up to 10"^ sites and a large number 
of different disorder configurations. We found that the 
Coulomb gap persists in the entire insulating phase but 
becomes narrower when approaching the MIT. The de- 
pletion of the DOS in turn leads to an enhancement of 
localization close to the Fermi level. This enhancement 
seems to be in contradiction to the results of Talamantes 
ct al. who inferred an interaction-induced delocaliza- 
tion from an investigation of the many-particle spectrum. 
Since the HF method is a severe approximation it is not 
clear whether the seeming disagreement is an artificial 
result of the HF approximation or whether it can be at- 
tributed to the use of different criteria Q for localization 
in a many-particle system or being in a different part of 
parameter space. 

In this paper we therefore investigate the range of va- 
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FIG. 1. Comparison of the exact and Hartree-Fock results 
for the single-particle DOS of a quantum Coulomb glass with 
3x4 sites for Wo = l, K — 0.5 and t — 0.1. The data represent 
averages over 1600 different disorder configurations. 



lidity of the HF approximation by comparing it to the 
results of numerically exact diagonalization of small lat- 
tices with up to 16 sites. 

We first discuss the single-particle (i.e. tunnehng) DOS 
g{e). In a many-particle system in general the DOS is 
defined via the single-particle Greens function 



g{e) = -- ImtrG^(£). 
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(2) 



Within the HF approximation the single-particle energies 
are simply given by the eigenvalues of the HF Hamilto- 
nian. A comparison of the single-particle DOS resulting 
from the HF and exact calculations for one particular set 
of parameters in Fig. ^ shows good agreement. We have 
performed analogous calculations for different values of 
the overlap {t = . . . 0.8) and filling factors {K = 0.5 and 
0.25). For all parameter sets investigated we found that 
the results of both methods agree within the statistical 
errors. Thus we conclude that the single-particle DOS of 
the quantum Coulomb glass is well described within the 
HF approximation. We note, however, that we cannot 
make a statement about the validity of the HF approx- 
imation asymptotically close to the Fermi energy since 
the small systems which we can study by exact diagonal- 
ization always possess a sizeable hard gap which obscures 
the Coulomb gap. 

We now turn to the localization properties. The sim- 
plest measure of localization for a single-particle state 
lip^) in a non-interacting system is the inverse participa- 
tion number 



(3) 



where the sum runs over all sites j. In practice it is often 
averaged over all states with a certain energy e 
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FIG. 2. Comparison of the Hartree-Fock and exact results 
for the return probability Rp for a particle in a quantum 
Coulomb glass with 3x4 sites for Wo = l, K = 0.5 and t = 0.3. 
The data represent averages over 1600 different disorder con- 
figurations. 
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The generalization of the participation number to many- 
particle systems is not straightforward, for a discussion 



see e.g. Ref. 



A consistent generalization of P is 



the probability Rp for a particle to return to its starting 
site in infinite time, which can be expressed in terms of 
single-particle Greens functions 

i 

For non-interacting electrons P~^{e) — Rp{e). In Fig. ^ 
we compare the return probability obtained by Hartree- 
Fock and exact calculations. Close to the Fermi energy 
the HF and exact results agree reasonably well. Away 
from the Fermi energy the return probability obtained 
by exact diagonalization decreases drastically which is 
not correctly described by HF jsj. We have carried out 
analogous calculations for different values of overlap (t = 
. . .0.8) and filling factors {K = 0.5 and 0.25). For all 
parameter sets we found the same behavior, viz. the HF 
method overestimates the return probability away from 
the Fermi energy whereas both methods agree close to 
the Fermi level. 

To summarize, we have compared the results of a HF 
approximation and of an exact diagonalization of the 
quantum Coulomb glass. We found that the results of 
the two methods for the DOS agree reasonably well for 
all energies. The localization properties are correctly de- 
scribed by HF close to the Fermi energy only. In par- 
ticular, the interaction-induced enhancement of localiza- 
tion which we found Q within the HF approximation is 
also present in the exact calculation. Consequently, it 
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appears that the seeming contradiction between the HF 
calculation and Ref. is not an artefact of the HF 
approximation but either due to the use of different quan- 
tities to characterize localization in a many-body system 
or due to being in different parts of parameter space. 
Thus a detailed investigation of the different localization 
criteria seems to be necessary The results of the 
present paper not only justify the application of the HF 
method to analyze at least some properties of the quan- 
tum Coulomb glass (taking into account the fact that the 
relevant physics at low temperatures is dominated by ex- 
citations close to the Fermi energy), they also suggest an 
improved calculational scheme. It consists of two steps: 
(i) solving the model within HF approximation and (ii) 
expanding the many-particle states in an energetically 
cutoff HF basis and diagonalizing this reduced problem. 
Such a scheme will enable us to treat comparatively large 
lattices and still obtain almost exact results for energies 
close to the Fermi energy. Work along these lines is in 
progress. 
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